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Abstract. We introduce a Markov Chain Monte Carlo simulation and data analysis 
package that extends the Cmbeasy software. We have taken special care in 
implementing an adaptive step algorithm for the Markov Chain Monte Carlo in order to 
improve convergence. Data analysis routines are provided which allow to test models of 
the Universe against measurements of the cosmic microwave background, supernovae la 
and large scale structure. We present constraints on cosmological parameters derived 
from these measurements for a ACDM cosmology and discuss the impact of the different 
observational data sets on the parameters. The package is publicly available as part 
of the Cmbeasy software at www.cmbeasy.org. 
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1. Introduction 

The wealth of recent precision measurements in cosmology [HElinilllEllSllZllHlini HIH 
ITT| 1121 Cni UH US] places stringent constraints on any model of the Universe. Typically, 
such a model is given in terms of a number of cosmological parameters. Numerical tools, 
such as Cmbfast [in], Camb jHj and Cmbeasy [TH], permit to calculate the prediction 
of a given model for the observational data. While these tools are comparatively fast, 
scanning the parameter space for the most likely model and confidence regions can 
become a matter of time and computing power. 

The cost of evaluating models on a n-dimensional grid in parameter space 
increases exponentially with the number of parameters. In contrast, the Markov Chain 
Monte Carlo (MCMC) method scales roughly linearly with the number of parameters 

EDI 1^ ■ The MCMC method has already been used to constrain various models 
121 121 123 123 • A popular tool for setting up MCMC simulations is the COSMO- 
Mc package |2I] for the Camb code, an improved proposal distribution for the local 
Metropolis algorithm has been proposed in 

In this paper, we introduce the AnalyzeThis package^ for Cmbeasy. It includes 
a parallel MCMC driver, as well as routines to calculate the likelihood of a model with 
respect to various data sets. We took special care in designing a step-proposal strategy 
that leads to fast convergence and mixing of the chains. This strategy is applied during 
the early stages of the simulation. As soon as the likelihood contour has been roughly 
explored, the adaptive step proposal freezes in. This ensures that the MCMC results 
are not contaminated by the adaptive steps. At the same time, the automated step 
optimization considerably improves performance and is rather convenient. The raw 
data files can be processed from within a graphical user interface (gui). Using the gui, 
one can marginalize, visualize and print one and two dimensional likelihood surfaces 
(see figure UJ). 

The plan of this paper is as follows: we describe the MCMC method and our 
implementation in section |2l A brief introduction to the software is given in section |H1 
We analyze the constraints on cosmological parameters from observational data sets in 
sectional In section El we present our conclusions, while the format of the MCMC data 
files is defined in [Appendix A 

2. Markov Chain Monte Carlo simulation 

In the following, we will review the basic ideas of Markov Chain Monte Carlo simulation 
[23 12H1 121] and describe our implementation. 

Suppose that 9 denotes a vector of all model parameters, X represents some 
observed data and L{X\6) is the likelihood of observing X given parameters 6. 
Specifying a prior distribution P{6) for the parameters, Bayes' theorem yields the 

J It is part of the cmbeasy v2.0 release. 
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Figure 1. The graphical user interface of Cmbeasy. It can be used to marginahze, 
visuahze and print the one and two dimensional likelihoods from the MCMC chains. 
Shown is the marginalized likelihood in the flmh^ ~ rif,/i^ plane of a ACDM model. 



posterior distribution tt{6\X) of 6 given the observed data X: 

= J Pie)Lix\e)de- 

Using 7t{6\X), one can compute expectation values 

as well as confidence levels. 

The idea of the MCMC method is to directly draw samples from the posterior 
7r(6'|X). The statistical properties of 7r(6'|X) may then be estimated using this sample. § 
To accomplish the sampling of 7r{6\X) one uses a Markov Chain, which is a stochastic 
process {9o,9i, ...,9n} where 9n only depends on 6'„_i. The idea is to choose the next 
point in the chain based on the previous point such that 7r(6'|X) becomes the stationary 
distribution of the chain 

Dist{^o, •••,^n} ^ vr(e|X) as n ^ oo. (4) 

There are several methods to accomplish this. We will concentrate on the Metropolis 
algorithm and its implementation in Cmbeasy. 

§ For instance, one can infer the mean of parameter 0*^'-' from the sample of parameter values 9i, ...TOn 
(the ergodic average): 

1 " 

= (3) 
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Figure 2. Illustrating the Metropolis algorithm for two parameters. Filled circles 
represent points belonging to the chain, empty circles are proposed but rejected points 
not belonging to the chain. In this example, the chain would be [9o, Oi, di, 63, 6*4, . . .]. 

2.1. The Metropolis Algorithm 

The algorithm is defined as follows (for an illustration see Fig. |2I): 

(i) Choose starting parameter vector 6*0. 

(ii) Compute the likelihood Lq{X\9o) of observing the experimental data given the 
parameters 6*0- 

(iii) Obtain a new parameter vector by sampling from a "proposal distribution" 
q{9i-i,9i) (see section EISl) • 

(iv) Compute the likelihood Li{X\9i). 

(v) If Li > Lj_i then save 9i as new point in the chain ("take the step") and go to (Imj). 

(vi) If Li < Li_i then generate a random variable u from [0, 1]. If m < Li/Li_i take the 
step as in ((^. If m > Li/Li_i then reject 9i, save 6'j_i as new point in the chain 
and go to 

This algorithm assumes fiat priors P{9) and a symmetric proposal distribution 
q{9i-i,9i) = q{9i,9i-i). Note that we assign likelihood zero to any parameter set that 
has at least one point outside its prior. In this version of the Metropolis algorithm, all 
parameters change with every step, a strategy called global Metropolis.! 

2.2. Convergence Testing 

At the beginning, the chain migrates from its random starting point to regions of higher 
likelihood. Points during this "burn-in" do not constitute a sample from 7r(6'|X) and 

f One may also change one parameter (or a subset of parameters) at a time. Such a local Metropolis 
algorithm is implemented in CoSMO-Mc|21). 
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should be eliminated. In principle, it may be difficult to tell from a single chain if it has 
converged towards the underlying 7r(0|X). In MCMC, one therefore uses several chains 
with random starting points and monitors mixing and convergence. 

The convergence test of Gelman and Rubin monitors the variance of a 
parameter between the chains. To be precise, consider using the last n points of each 
of m chains for the test. Let ipij label one entry of the parameter vector 6 at point 
j = 1, . . . , n in chain i with ip^ the mean for chain i and ip the mean of all chains. The 
variance between chains B and the within-chain variance W are then given by 

m 

1=1 

^ = - E -here s^, = -— ^ ^^(V^., - (6) 

i=l j=l 

and the quantity 



R = ^^^r^ (7) 

should converge to o\ie.\ A value of i? < 1.2 for all parameters indicates that the chain 
is sampling from '7r(^|X).§ From this point onwards one may use the chain points for 
parameter estimates. 

When do we have enough points for parameter estimation? This question is not easy 
to answer, since it is depends on the model, the used data sets and the desired accuracy 
how many chain points are needed for a robust estimate of parameters. Therefore, in 
our implementation the MCMC simulation runs indefinetly. However, any "breaking- 
criterion" may be implemented easily, and the chains may be monitored with external 
programs during the run. 



2.3. Adaptive Step Size Gaussian Sampler 

The number of steps needed for good convergence and mixing depends strongly on 
the step proposal distribution. If the proposed steps are too large, the algorithm will 
frequently reject steps, giving slow convergence of the chain. If, on the other hand, 
the proposed steps are too small, it will take a long time for the chain to explore the 
likelihood surface, resulting in slow mixing. In the ideal case the proposal distribution 
should be as close to the posterior distribution 7r(6'|X) as possible - which unfortunately 
is not known a priori. While a simple Gaussian proposal distribution with step sizes 
CTfc is sufficient, it is not optimal in terms of computing costs if cosmological parameters 
are degenerate. A naive Gaussian sampler would move slowly along the degeneracy 
direction, unaware of any degeneracy (see figure 0)). 

I In a realistic situation, the numerator is an overestimate whereas the denominator is an underestimate 
of the variance of the stationary distribution of %p. 

§ Indicating that the step sizes and directions wiU freeze in in our implementation. 
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Figure 3. Illustrating the naive Gaussian sampler with fixed step sizes for two 
parameters. The (unknown) true likelihood surface is shown in black, the proposal 
distribution with arrows in red. This proposal distribution does not does not take into 
account the degeneracy among the parameters and 0'^', leading to slow mixing. 



Instead of using a naive Gaussian proposal distribution, we sample from a 
multivariate Gaussian distribution with covariance matrix estimated from the previous 
points in the chains. By taking into account the covariances among the parameters, 
we effectively approximate the likelihood contour in extent and orientation - the 
Gaussian samples are taken along the principal axis of the likelihood contour. Denoting 
6i^i — 6i = Ui for notational convenience, the proposal distribution we use for the steps 
is 

1 



~ ATexp 



-U' S Uj 



(8) 



(9) 



where N = (27r) ' (det S) and S is the covariance matrix 
( o-l pi2 ... PlK \ 

2 

P2I <^2 ■ ■ ■ P2K 

\ PKI ■ ■ ■ Pk-ik (^"k / 
The sampling is most easily performed by diagonalizing the covariance matrix 

T = D T^S-^T = D-\ (10) 

where T is an orthogonal matrix. Using this, equation (jH)) becomes 

g(^,_i,^,) ~ iVexp [-ujTT'^S-^TT^Uil (11) 



A^exp 



-v! D 



Vi 



(12) 



where Thus, the procedure for obtaining a sample Ui is as follows: 
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(i) Find the eigenvalues (t| and eigenvectors of S. Construct the transformation matrix 
T from the eigenvectors. 

(ii) Draw Gaussian samples with variances aj, thereby obtaining the vector Vi. 

(iii) Then itj = Tvi is the desired sample from the multivariate Gaussian with covariance 
matrix S. 

The convergence can be further improved by scaling the covariance matrix S with 
a variable factor a. Using a, we can cope better in situations where the projected 
likelihood takes on banana shapes such as in |21]. It also improves the convergence 
during the early stages when the low number of points available limits the estimate of 
the covariance matrix S. We dynamically increase a if a chain takes steps too often, | 
while we decrease a if the acceptance rate is too low.§ The positive effect of our scheme 
on the convergence is illustrated in figure |3] 

One can show that modifying the proposal distribution based on previous chain data 
during the run may lead to a wrong stationary distribution 7c'{d\X) fI^ \'2S\. Therefore, 
we only apply the dynamical strategy of finding an optimal step proposal during the early 
stages of the simulation. When the convergence is better than R = 1.2 and the chain 
has calculated a certain number of points, we freeze in the step proposal distribution. 
All points in the chain before this freeze in should be discarded (see also [Appendix A ). 

3. The Software 

The package is part of Cmbeasy and consists of two main components. The first one is 
a MCMC driver using LAM/MPI [32j for parallel execution of each chain. The second 
one is the AnalyzeThis class which is designed to evaluate the likelihood of a given 
model with respect to various data sets. These sets include the latest data of WMAP 
TT and TE jSlE!, ACBAR CBI 0, VSA (Hj , 2dFGRS [HI 1121 CHI, SDSS HHIini, 
the SNe la compilations of Riess et al. [7], Tonry et al. |H] and Knop et. al. [H] as well 
as the IfA Deep Survey SNe la data ^0] • Data files for all experiments are included for 
convenience. New data sets are added continuously to the code. 

3.1. The MCMC Driver 

The example MCMC driver consists of two routines: master() and slave(). Using 
LAM/MPI for parallel computing, one master and up to ten slaves may be started. 
The master() will determine the initial random starting position for each chain. In a 
never ending loop, it then sends the parameters to the slave()'s and collects the results 
when the computation is finished. Whenever a step has been successful, it stores the 
parameters and likelihoods of the last step together with the number of times the chain 

I Frequent acceptance means that the hkehhood at the next step is roughly comparable to the current 
one. This happens when the chain rarely takes steps larger than la. 

§ Rare acceptance means that the chain rarely explores points with the same likelihood, i.e. neighboring 
points within la. 
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Figure 4. Convergence properties using different proposal distributions for a 
cosmological model with seven parameters. For illustration purposes we display the 
average i?-statistic as a function of number of models computed. The univariate 
Gaussian proposal distribution using fixed variances but adaptive overall step size a 
(black dashed curve) shows convergence after about 2800 iterations. The multivariate 
Gaussian proposal distribution with covariance matrix estimated from the previous 
chain points and adaptive overall step size a as suggested in this paper (red solid 
curve) converges after 500 iterations. 



remained at the same point in parameter space in a new line of one file per chain. The 
format of these raw MCMC data files is defined in [Appendix A[ The master () monitors 
convergence and mixing and determines the next step for the slave(). Before freeze in, 
the covariance is estimated and the step proposal accordingly modified. After freeze in, 
the proposal distribution remains unchanged. 

3.2. The AnalyzeThis Class 

The AnalyzeThis class provides several routines concerning CMB, SNe la and Large 
Scale Structure measurements. It also contains routines for marginalizing and plotting 
the Monte Carlo data. 



3.2.1. WMAP The WMAP routines are (slightly modified) C++ ports of the 
likelihood code ^ available at the LAMBDA ^ web site. When the WMAP 
routines are called for the first time, the covariance matrices provided by WMAP 
will be converted to a binary format to speed up future use. A routine for WMAP 
normalization of the Ci spectrum using the binned TT may be used instead of the old 
COBE normalization. (After the quick normalization, a best fit normalization may be 
called which uses the full likelihood routine provided by .25 .) 
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3.2.2. ACBAR, CBI and VSA We use the procedures for likelihood computation 
described by the ACBAR, CBI and VSA collaborations jHEUHl, using window functions 
and calibration uncertainty. One can de-select data bins from each of these datasets. It 
is, for instance, possible to calculate the likelihood for the low / data only, thus speeding 
up the computation of a model, because the window functions of the high / data require 
multipoles up to Z = 4000. One may also wish to exclude I < 800 data if one is using 
the WMAP observations in order to keep the sets independent. 

3.2.3. 2dFGRS For this Large Scale Structure dataset one compares the data with the 
theoretical power spectrum at z = 0.17, the effective redshift of the survey, multiplied 
with the window function. We only include the region with k/h < 0.15 Mpc~^ since at 
smaller scales nonlinear effects need to be taken into account jUj. For these values, the 
bias is nearly constant ^21 El- One may either specify the bias or marginalize over it. 

3.2.4. SDSS The theoretical power spectrum to be compared with the data should 
be evaluated at the effective redshift of this survey, at z = 0.1. We use the data 
given in Table 3 of J3] and the appropriate window functions. One may select the 
maximum fc-value to be included in the likelihood estimate, but including data beyond 
k/h > 0.15 Mpc~^ requires nonlinear corrections. Again, one can specify the bias or 
marginalize over it. 

3.2.5. Supernovae la We include four routines for calculating the likelihood with 
respect to SNe la data. Please note that the sets of Riess et al., Tonry et al. and 
Knop et al. are not independent. 

Riess et. al. One can use the full dataset, subset selection of the "gold" set as described 
in [Z| is possible. Likelihood computation as given in this paper. 

Tonry et. al: From the supernovae compilation of Tonry et al. [S] one can use the 
full data set of 230 Supernovae or, alternatively, one may use a restricted set of 172 
supernovae, where supernovae with z < 0.01 and with excess reddening have been 
omitted as suggested in |H|. In any case, we have taken the particular velocity uncertainty 
to be v = 500 km/s corresponding to Az = 0.00167 and computed the corresponding 
uncertainty in the luminosity distance to obtain the likelihood: 



Here, ci^^^ and rfp is the experimental and theoretical luminosity distance, respectively. 




N 



(13) 



(14) 
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Table 1. Flat priors on the parameters used in our MCMC simulations 



Parameter Min Max 





0.016 


0.03 




0.05 


0.3 


h 


0.60 


0.85 


T 





0.9 


ris 


0.8 


1.2 



Knop et. al: The 54 SNe la presented in Tables 3-5 have been included in the code. 
We use the data with strech and extinction correction applied, subsample selection as 
discussed in P is also possible. 

If ADS survey: One may choose to include all 23 SNe la given in ^U], or one may 
exclude the supernovae with excess reddening and those not unambigously identified as 
SNe la as described in [TU]. 

3.3. The Graphical User Interface 

The gui may be used to process the raw output files of the MCMC chains. After 
starting cmbeasy, the first step is to "distill" the chain data files. By distilling we 
mean the merging of the raw chain output files into one file. In addition, distilling 
removes all burn-in data. To get started immediately, we include raw data from an 
example MCMC run in the resources directory of cmbeasy. The four chains are called 
montecarlo-chaina.dat. They are runs for a ACDM model and need to be distilled first. 
Two and one dimensional marginalized likelihoods may then be plotted and printed from 
within the gui (see figure Q). Please see the "howto-montecarlo" document shipped with 
cmbeasy for an introduction. 

4. Cosmological constraints from the data sets 

Having described the method and software, we can now proceed to investigate the impact 
of the different data sets on the distribution of parameters for a given cosmological 
model. For illustrative purposes we limit ourselves to a flat ACDM cosmology with 
flve parameters. We take the reduced baryon and matter densities flbh"^ and VLmh"^, the 
Hubble parameter h = ifo/(100km s~^Mpc~^), the optical depth to the last scattering 
surface r and the spectral index of the initial power spectrum Ug as parameters. We 
neglect any tensor contributions and we marginalize over the amplitude of the initial 
power spectrum. Thus, the amplitude is treated as a "nuisance" parameter that is 
integrated out. 

In table [T] we display the cosmological parameters and the flat priors used. The 
limits of the parameter ranges are chosen in accordance with previous results Pl E] • 
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Figure 5. CMB data set for constraining cosmological parameters. 



First, we will investigate constraints from CMB-only data, subsequently adding 
Supernovae la and large scale structure data sets. For a detailed description of how 
we implemented the datasets we refer the reader to section 13.21 of this paper. We have 
performed three MCMC simulations with CMB only, CMB+ SNe la and CMB + SNe 
la + LSS data sets. Convergence was reached in each run after about ~ 2000 iterations, 
each simulation was run until 55 000 models had been computed. 

For the CMB-only run we use the measurements of WMAP (TT and TE spectra), 
CBI, VSA and ACBAR up to / = 2000, removing data points from CBI, VSA and 
ACBAR where we have WMAP measurements with comparable error bars to keep the 
data sets independent. The entire set we used is displayed in figure El For CMB + 
SNe la we add the data set of Riess et al. ("gold sample"), and finally for CMB + 
SNe la + LSS we add the large scale structure measurements of the SDSS collaboration 
using all points with k/h < 0.15 Mpc~^ where the perturbations are still linear |15j . 
The results are presented in table El marginalized likelihoods for the parameters are 
displayed in table 01 and figure IHl 

We will now discuss the impact of the different sets on the parameters. Consider 
figure ini The CMB-only data set does not well constrain Qmh"^- Adding supernovae 
improves the bounds on the total matter content considerably, because the lumiosity 
distance depends sensitively on Qm- Adding large scale structure data breaks the 
degeneracy between matter contribution and initial power spectrum amplitude, hence 
we obtain still tighter bounds. Note that the distribution of Qmh"^ obtained from CMB- 
only and CMB + SNe la -f LSS data sets differ somewhat. This has already been noted 
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Figure 6. One dimensional marginalized distributions for the parameters. Constraints 
from CMB only (red, dashed), CMB + SNe la (blue, dotted) and CMB + SNe la + LSS 
(black, solid). The histograms extracted from the chains were fitted using a function 
f{x) = exp(p(a;)) with p{x) a sixth-order polynomial jl4| . 



Table 2. Constraints on the parameters from a combination of data sets. 
These confidence intervals were generated from the one-dimensional marginalized 
distributions (errors are given at 68.3 % confidence level). 



Parameter CMB only CMB + SNe la CMB+ SNe la + LSS 





0.0226liro^}? 


U.UZZO<_Q QQPg 


U-U^^'O-0.0009 




0.1306t°:l]}^ 


u.ioao_Q QQg 




h 


0.7161^04? 


u.oa»_g 


0.6821°:°- 


T 


0.1019iSJ:°f 


0941+° °^® 


U.U900_Q QQ43 


ris 




n Qc;'?+0-024 


U.900_o.019 



in P3] . Even from CMB measurements alone we can infer the presence of dark matter at 
high significance. In contrast, the baryon contribution Qhh'^ is already well constrained 
by the CMB alone, adding SNe la and LSS data only improves the bounds slightly. 
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Table 3. Two-dimensional marginalized likelihoods for the cosmological parameters 
using CMB + SNe la + LSS data. The contours are one, two and three sigma 
confidence regions, respectively. 




0.13 0.14 0.15 0.16 




0.13 0.14 0.15 0.16 




0.13 0.14 0.15 0.16 




0.13 0.14 0.15 0.16 



0.026 
0.024 
^ 0.022 



0.02 




0.62 0.66 0.7 0.74 0.78 

h 





0.62 0.66 0.7 0.74 0.7S 

h 





0.62 0.66 07 0.74 0.78 

h 



0.9 0.93 0.96 0.99 1.02 1.05 



The CMB is most sensitive to this parameter since what one observes are essentially 
oscillations in a photon-baryon plasma, the density of baryons is a critical parameter 
for the shape of the CMB spectrum. SNe la and LSS measurements are more or less 
insensitive to this parameter, though hints at oscillations in the power spectrum have 
been detected [TT] . 

The Hubble parameter h, is only slightly constrained by CMB measurements alone. 
Adding SNe la data considerably improves the bounds, adding large scale structure data 
does not improve the bounds very much. The bounds obtained for the Hubble parameter 
for the full data set are consistent with the HST key project value h = 0.72 ± 0.08 
derived from measuring the local Hubble flow. Even though the optical depth does not 
directly influence SNe la predictions, adding SNe la data tightens the bound on the 
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optical depth r as seen in figure IHl The mechanism is somewhat indirect: SNe la data 
tightens the bound on Q.^ and thus hmits the range of allowed values for r, h and n<, 
from CMB measurements. As we marginalize over the bias of the SDSS data, large scale 
structure does not add further to the bounds on r. 

5. Conclusions 

We have introduced the AnalyzeThis package, which can be used to constrain 
cosmological models using observational data sets. The AnalyzeThis class provides 
many functions to compute the likelihood of a model with respect to measurements of 
the CMB, SNe la and Large Scale Structure. jj 

In order to constrain models of the Universe with a substantial number of 
parameters, we include a Markov Chain Monte Carlo driver. As the MCMC step 
strategy determines the convergence speed of the chains, we implemented a multivariate 
Gaussian sampler with an additional dynamical scaling. We stop the adaptive 
improvement of the step proposal density as soon as a good level of convergence is 
reached and discard all data calculated before. Using this approach, we combine the 
advantage of a single-run automatic optimization of the step strategy with the demand 
of a static step proposal density. 

The output of the Monte Carlo chains is in a human readable format and may be 
processed by any software, even during the run. For convenience however, one may 
use Cmbeasy's gui to marginalize, plot and print two and one dimensional likelihood 
surfaces. 

Finally, we discussed the impact of the different data sets on the parameters for 
the case of a ACDM cosmology with five parameters. 

Acknowledgments 

We would like to thank Robert R. Caldwell, Robert A. Knop, Havard B. Sandvik, Max 
Tegmark, Adam G. Riess and Timothy J. Pearson for helpful discussions. M. Doran 
is supported by NSF grant PHY-0099543. CM. Miiller is supported by GRK grant 
216/3-02. 

Appendix A. Files and Formats 

Appendix A.l. MCMC chain output data 

The parameters and likelihoods of each model are stored in files with the naming con- 
vention "montecarlo_chaina.dat", where "a" is an integer. Each file corresponds to a 
slave() process (and hence to an independent chain). Each line in such a file represents 

tt We would be happy to include contributions of routines for recent and future measurements. So if 
the readers favorite measurement is not included yet, please send in a few lines of code. 



Analyze This! A cosmological constraint package for cmbeasy 
a new point in parameter space. The format of one such hne is: 



15 



PO \ Pi \ ■ ■ -IPn \ L \ lo \ . ..\ lm\ ai \ . . . I Ofc I M, 

where Pi are the parameters of the model, L is the overall likelihood, the /j are likelihoods 
of different experiments, are auxiliary fields (to store as, for instance). Finally, M 
determines the weight of this parameter set, i.e. the "time" spend before leaving this 
point in parameter space. When the step sizes and the covariance estimates are frozen in, 
one line in each file gets a zero M. All data before this line should be regarded as "burn- 
in". The graphical user interface (using the distillChains() routine of AnalyzeThis) will 
automatically discard all models before freeze-in. 

Appendix A. 2. Monitoring Progress and Covariance 

The master() routine outputs the convergence and some more information into a file 
called "progress.txt". Until the step proposal distribution is frozen in, the covariance 
matrix is in addition output to "covMatrix.txt". Messages of any errors occuring in 
slave() will be appended to the file "errorlog.txt". 
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